Complete homochirality induced by the nonlinear autocatalysis and recycling 
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A nonlinear autocatalysis of a chiral substance is shown to achieve homochirality in a closed 
system, if the back-reaction is included. Asymmetry in the concentration of two enantiomers or 
the enantiometric excess increases due to the nonlinear autocatalysis. Furthermore, when the back- 
reaction is taken into account, the reactant supplied by the decomposition of the enantiomers is 
recycled to produce more and more the dominant one, and eventually the homochirality is estab- 
lished. 

PACS numbers: 



Natural organic molecules associated with living mat- 
ters usually have two possibilities in their stereo- 
structures, in a right-handed form(R) or in a mirror- 
image left-handed form(S)i£ These two forms are called 
enantiomers, and the molecules are said to be chiral. Chi- 
ral molecules are optically active to rotate the direction 
of polarization of plane polarized light. From the en- 
ergetic point of view, these two enantiomers can exist 
with an equal probability, but the life on earth utilizes 
only one type: only levorotatory(L)-amino acids (or S) 
and dextrorotatory (D)-sugards (or R). This symmetry 
breaking in the chirality is called the homochirality. 

The origin of this unique chirality has long intrigued 
many scientists^ In order to find the physical origin 
of this homochirality, initial asymmetry in the primor- 
dial molecular environment has to be created by chance 
or engendered deteminately by external or internal fac- 
tors, such as the parity breaking effect in the weak 
interaction^^& the asymmetry in circularly polarized 
lightj2i£ or adsorption on optically active crystals^ Then 
the induced small initial chiral asymmetry has to be am- 
plified. 

Frank has shown theoretically that an autocatalytic 
reaction of a chemical substance with an antagonistic 
process can lead to an amplification of enantiometric ex- 
cess (ee) and to homochiralityiS Many theoretical mod- 
els are proposed afterwards, but they are often criti- 
cized as lacking any experimental support i Recently, 
asymmetric autocatalysis^ of pyrimidyl alkanol has been 
studied intensivelyiii*i2iiii4 The enhancement of ee was 
confirmed^ and its temporal evolution was explained 
by the second-order autocatalytic reactionii&M But only 
with the nonlinear autocatalysis, chirality selection is not 
complete and the value of ee stays less than unity. Here 
we show that the complete homochirality is achieved by 
the back-reaction to recycle the reactant. If the rate of 
back reaction is small, it takes a long time before the 
homochirality is achieved. 

We consider a chemical reaction such that substances 
A and B react to form substance C. Though reactants A 
and B arc achiral, the product C happened to be chiral in 
two enantiometric forms; i?-isomer (i?)-C and S'-isomcr 
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In the experiment by Soai et al4^ one provides a small 
amount of enantiomers with concentrations Rq and So re- 
spectively in the reactants of concentrations A$ and Bo- 
Instead of the open flow system considered by Frank^ 
we assume a closed conserve system such that concentra- 
tions R and S at time t vary in proportion to the present 
amount of the reactants A and B as 

^ = k R (A + Ro + Sa-R-S)(B + R Q + S -R-S) 
at 

dS 

— = k s (A + R0 + S0-R- S)(B + R a + S a -R-S). 

(2) 

Here we first neglect the back reaction from (S)-C or 
(R)-C to A and B. Later, the effect of back reaction is 
shown to be very important in achieving the complete 
homochirality. 

The positive reaction coefficients, kn.s, may depend 
on the concentrations R and S, if the reaction is auto- 
catalytic. We consider various cases step by step in the 
following. Furthermore, for a simplicity reason we as- 
sume that there is an ample amount of substance B, and 
the reaction is controlled by the less abundant substance 
A. Then by normalizing the concentration by the initial 
concentration A + R + So as r = R/{A + R + So), 
s = S/(Ao + Rq + Sq), the equation is simplified as 



dr 

~dt 
ds 

~dt 



= k r (l — r — s) 
k s (l — r — s) 



(3) 



where the reaction rate coefficient k r _ s = k^s(Bo + Ro + 
So) may depend on r and s. 

Since both enantiomers arc equivalent except the hand- 
edness, the time evolution should be symmetric in r and 
s. The chirality of the system can be brought about 
through the asymmetry of the initial concentrations ro 
and So of two enantiomers, (R)-C and (S)-C. In the fol- 
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lowing wc study the temporal variation of the enatiomct- 
ric excess measured by 



n = 



(4) 



Non-Autocatalytic Case 

If the reaction is not autocatalytic, the reaction coef- 
ficients are simply constant as k r = k s — k. Then the 
ratio of the two equations (3) gives dr/ds = 1, which de- 
fines the trajectory in the r — s space. It is easily solved 
as r — s = ro — Sq =const. The time evolution of enan- 
tiomers is depicted in the flow diagram in Fig. 1(a). The 
flow is on the line with a unit slope, and terminates on a 
point of crossing the diagonal line r + s = 1; the reaction 
stops because the reactant A is consumed up in a closed 
system. The diagonal line r + s = 1 is in fact a set of 
stable fixed points. 
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(c) Nonlinearly autocatalytic (d) with back reaction 

FIG. 1: Flow diagrams in r-s space, (a) Nonautocatalytic, 
(b) linearly autocatalytic, (c) nonlinearly autocatalytic and 
(d) nonlinearly autocatalytic with backward reactions (A/fc = 
0.1). 

The ee decreases since the numerator r — s in fi re- 
mains constant whereas the denominator r + s increases 
in time. Actually, the equation is easily solved, and ee is 
obtained as O = (r - s )/(l - (1 - r - s )e~ 2kt ). The 
ee decreases from the initial value fio = ( r o — sq) I ( r o + s o) 
to the asymptotic value fioo = r$ — sq, which should be 
very small if the initial amounts ro = Ro/ (Aq + Rq + Sq) 
and so = Sq/(Aq + Rq + So) are small. This corresponds 
to racemic mixtures. 

Linearly Autocatalytic Case 

If the reaction is autocatalytic in the first order, the 
reaction coefficient is linearly proportional to the product 
concentrations as k r = kr and k s — ks. In this case, the 
ratio of two rate equations gives dr/ds = r/s, which is 



easily solved as r/s = ro/so =const. The time evolution 
is depicted in the flow diagram Fig. 1(b). The flow is 
radial from the origin and trivially the ee does not change 
even as the reaction proceeds, 



fi 



^0 ~ Sq 

ro + s 



fin 



(5) 



as first notified by Franks 

Nonlinearly Autocatalytic Case 

If the reaction is autocatalytic in the higher order, the 
reaction coefficient depends nonlinearly on the product; 
for the second order autocatalysis, for example, k r = kr 2 
and k s = ks 2 . In this case, the ratio of two rate equations 



gives 



dr 
ds 



(6) 



which is easily solved to give the trajectory in the r — s 
space as 1/r — 1/s = 1/rn — 1/sn =const. The time 
evolution is depicted in the flow diagram Fig. 1(c). Below 
the diagonal r = s, the flow bends down, whereas above 
the diagonal the flow bends up. The final concentrations 
Too and Soo should satisfy the relation roo + Soo = 1. Then 
the final value of the ee is determined by the initial value 
fin and the amount of the input ro (or sq) as 



fic 



sgn(fio) 




1-fin 



nr 



Ifin 



(7) 



If one starts with very small values of ro and So, the ee 
approaches | fioo | « 1— 7*o(l— fio)/|fio|- The amplification 
of ee is expected to be large if the initial composition ro 
and so are small. For example, with the same initial ee as 
fi = 2%, if one starts with r = 0.051 and s = 0.049, 
the final ee is fioo = 19% according to eq.(6), whereas 
with the initial values r = 0.00051 and s = 0.00049, 
then fioo = 98%. 

Homochirality by Back Reaction 

From the previous analysis we found that the ee am- 
plification takes place with a nonlinear or higher-order 
autocatalytic chemical reaction. However, the final chi- 
rality is not complete. The ee is smaller than 100%, even 
the starting values ro and so are very small. Here we 
show that the inclusion of the back reaction from the 
products (R)-C and (S)-C to A brings about the drastic 
change in the flow diagram. 

A simple back reaction process is incorporated into the 
rate equation (3) as 



dr 

~di 
ds 



kr 2 (l — r — s) — Xr 
= ks 2 (l -r-s)- As 



(8) 



with a rate constant A of the back reaction. Now the 
flow diagram in r — s space changes from Fig. 1(c) to (d). 
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The diagonal line r + s = 1 is no more a set of infinitely 
many stable fixed points. Instead, there are seven fixed 
points in the flow diagram, Fig. 1(d). The trivial fixed 
point at the origin 0:(r o , s ) = (0, 0) is stable. On the r- 
axis there are two fixed points; Ui:(l — yl 4A/fc)/2, 0) 
is unstable and + - 4A//c)/2,0) is stable. There 
are analogous points on the s-axis as the unstable 
fixed point U 2 :(0,1 - y/l - 4A/fc)/2) and the stable 
fixed point S 2 :(0, 1 + y/l - 4A/fc)/2). On the diago- 
nal line, r = s, there are two other unstable fixed 
points: U 3 :((l - ^\ - 8A/fc)/4, (I - y/l - 8A/fc)/4) and 
t/ 4 :((l + V 1 - 8A/fc)/4, (1 + -8A/fc)/4). 

If the rate of back reaction A is small, t/1.2,3 are close 
to the origin O, and U4 and Si,2 are close to the diagonal 
line r + s = 1. Then the flow diagram looks similar to 
Fig. 1(c). However, there is a big difference in the flow 
diagram Fig. 1(c) and (d). In Fig. 1(c) the diagonal line 
r + s = 1 is a set of infinitely many fixed points. On this 
line, neither r nor s varies in time. On the contrary, in 
Fig. 1(d) there are only three stable fixed points; O and 
Si ,2- The trivial fixed point O corresponds to the struc- 
tureless or lifeless world, without the chemicals (R)-C or 
(S')-C. Its basin of attraction is a small region enclosed 
by curves connecting O and Ui^s- If the the priomor- 
dial environment allows some chemical system to reach 
exterior of this basin, chemical complexes with chirality 
can be produced. The nonlinear autocatalysis amplifies 
the ee, and afterwards the back reaction process selects 
only one type of a chiral component. 
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FIG. 2: Autocatalytic evolution of concentrations r (contin- 
uous curve) and s (broken curves) of two enantiomers with 
(A = 2 x 1CP 4 ) and without (A = 0) back reaction. 

The time evolution of the concentration of each com- 
ponent r and s is depicted in Fig. 2. The initial ee is 
as small as fi = 0.01% with r = 0.00050005 and 
s = 0.00049995. If there is no back reaction (A = 0), 
then in the early stage reactant concentrations increase 



as r « r o/(l — r ot) and s ~ So/(l — sgt), and tend 
to diverge at times l/r and 1/sq respectively. Before 
the divergence, the shortage of the reactant A leads to 
the saturation of r and s, as r rs ^(1 — e~ 2kt ) and 
s « Soo(l — e~ 2kt ). The final amounts and Soo are 
finite and nonzero. The ee is amplified to fioo =7.9% 
which is less than unity. 

If there is a small but a finite back reaction process 
with A = 0.0002, both r and s increase again and ap- 
proach close to the values corresponding to the unstable 
fixed point U4. Then the component which has a slight 
advantage, r in Fig. 2, starts to dominate, and the other 
chiral type extinguishes gradually. The substance (S)-C 
decomposes and supplies the reactant A which is con- 
summed to increase its enatiomer (i?)-C. The reactant A 
is, so to say, recycled. Eventually, the complete homochi- 
rality with fl =100% is achieved. 

We have shown theoretically that in a closed chemical 
system, the nonlinear autocatalysis amplifies the initial 
small enantiomctric excess. But eventually, the simple 
back reaction promotes the decomposition of less abun- 
dant enantiomer to the reactant, which is recycled to 
produce the more abundant type. Through this recycling 
process, the complete homochirality can be achieved. 

The ratio of the forward and the backward reaction 
rates, k and A, is related to the equilibrium concentra- 
tion or the yield. Therefore, a system with a yield less 
than unity might have sufficiently large back reaction rate 
to achieve complete homochirality in experimental time 
scales. But if the rate of back reaction is too large and 
the forward reaction is autocatalytic, the system might 
be trapped into the basin of attraction of the origin and 
the reaction can not proceed. On the contrary, when 
the rate of back reaction is too small, the enatiometric 
amplification can take place, but the final selection of chi- 
rality should take a long time. This might be the reason 
that the previous experiments by Soai et alii^ii are well 
fitted by a model with only the forward autocatalytic re- 
action. To observe the effect of back reaction, one might 
have to wait a long time, depending on the back-reaction 
rate. Otherwise one has to search appropriate chemical 
systems with proper reaction constants. 

If one considers spatial dependence of the chemical re- 
action, further complications arise in the process of chi- 
rality selection. Initial productions of enantiomers actu- 
ally take place randomly in space and time. The nonlin- 
ear autocatalysis promotes the chemical reaction around 
these initial enantiomers, and the reaction fronts propa- 
gate and expand. When fronts of different enantiomers 
meet with each other, the competition should take place 
between them. This competition of enantiomers in space 
is an interesting future research subject. 
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